Laser-induced quantum pumping in graphene 
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We investigate non-adiabatic electron pumping in graphene generated by laser irradiation with linear polariza- 
tion parallel or perpendicular to the transport direction. Transport is dominated by the spatially asymmetric 
excitation of electrons from evanescent into propagating modes. For a laser with parallel polarization, the 
pumping response exhibits a subharmonic resonant enhancement which directly probes the Fermi energy; no 
such enhancement occurs for perpendicular polarization. The resonance mechanism relies on the chirality of 
charge carriers in graphene. 
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The periodic modulation of a mesoscopic conductor 
attached to electronic reservoirs can be used to induce a 
transfer of electrons across the system, which produces 
a dc current in absence of a bias voltage. Such quantum 
pumps can be driven adiabatically slowly, provided that 
their work cycle is oriented, 1 which requires variation of 
at least two parameters. For faster, non-adiabatic driv- 
ing, single-parameter pumping is possible, 2 4 which has 
been observed in a number of experiments. 5 11 

With the advent of graphene, it has been suggested to 
exploit and probe the properties of this material both by 
adiabatic 12 16 and non-adiabatic pumping, 17 ' 18 with the 
driving typically provided by applying ac electrostatic 
voltages. In this work, by contrast, we consider driving 
graphene by linearly polarized laser irradiation, 19,20 as is 
depicted in Fig. 1. We show that this results in a sizable 
pumped current, which moreover is resonantly enhanced 
when the polarization of the electric field is parallel with 
the transport direction. The resonance occurs at a Fermi 
energy equaling half the photon energy, and thus may be 
useful to characterize the properties of a graphene sam- 
ple. The polarization dependence arises from a selection 
rule which directly probes the chirality of the charge car- 
riers in graphene. 

We consider a graphene monolayer sheet of length 2L 
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FIG. 1. Sketch of a graphene sheet of length 2L and width W 
connected to two highly doped graphene leads. The left half 
of the sheet is driven periodically in time by laser irradiation 
of frequency uj and amplitude Eq. 



and width W connected to two highly doped leads, as 
shown in Fig. 1. In the absence of the laser irradiation, 
the potential energy in the sheet is Ub(x) = 0, while in 
the leads Ub(x) <C 0; the potential energy is indepen- 
dent of the transverse coordinate y. The low-energy ex- 
citations are then described by the quasi-one-dimensional 
Dirac Hamiltonian Hq = hvpk • a + Ub(x), where 
vp ~ 10 6 m/s denotes the Fermi velocity of graphene. 
The components of a are the Pauli matrices <j x and a yi 
which act in a pseudo-spin space that physically describes 
the amplitudes on the two inequivalent lattice sites of 
graphene's honeycomb lattice. Because of the term k ■ cr 
in i^o, the electronic states of graphene possess chiral- 
ity, which is responsible for many of its properties, 21 
including those that result in efficient pumping by ac 
gating. 12 ' 17 

The pumping is induced by laser irradiation focused 
on a region of size L next to the left electrode, with lin- 
ear polarization along the unit vector n in the plane of 
the graphene sheet. The associated oscillating electric 
field E(x,t) = Eosm(ujt)a(x)n is modeled by a piece- 
wise constant laser spot profile a(x) (the magnetic field 
in the graphene plane has a negligible physical effect). 
To avoid complications arising from a position-dependent 
scalar potential, we work with the Weyl gauge in which 
the electric field is derived from a time-dependent vector 
potential, E(x,i) = —dA(x,t)/dt. Using the minimal 
coupling k — » k + eA for the electron charge — e, we 
obtain the time-dependent Dirac Hamiltonian 

H(t) = H + a(x) cos(ujt)U n , (1) 

where the constant coupling operator is given by 

Un = evp—n ■ a. (2) 

It follows that the polarization direction n couples to the 
pseudo-spin a and thus affects the quasiparticle chirality. 

The pumped current can be obtained by Floquet scat- 
tering theory, 22 which is valid as along as interactions 
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can be neglected. At low temperatures, the dc current is 
then given by the generalized Landauer formula 2 ' 23 



4e 



E F 



deAT(e), 



(3) 



where the factor 4 accounts for spin and valley degen- 
eracy. The central quantity is the net transmission im- 
balance AT(e) = £jr££(e) - T^(e)], which involves 

the transmission coefficients T^\e) from lead i — L,R 
to lead £' 7^ i, where n labels sidebands reached by ab- 
sorption or emission of n photons (n < corresponds to 
emission). The definition of the transmission coefficients 
includes the summation over the transverse wave num- 
ber k yi which here is conserved, and can be considered 
continuous for wide and short pumps (W ^> L). The 
pumped current is non-zero since the specified setup is 
asymmetric. In particular, as in the case of driving by an 
ac gate voltage, 12 ' 17 carriers entering the driving region 
from the left lead in evanescent modes can be promoted 
to propagating modes, which allows them to reach the 
right lead, while electrons arriving in evanescent modes 
from the right lead practically never reach the driven re- 
gion; the result is net transport from left to right. 

Before presenting quantitative results we establish the 
preferred pumping regime and the typical magnitude of 
the pumped current. There are two relevant parameters. 
The first is the effective driving strength 
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Adjacent sidebands possess relative energy shifts fkj, and 
are coupled by a "hopping" energy \U n \/2. The ratio of 
these two quantities, ^/p, equals the number of sidebands 
that get populated in the scattering process, 17 i.e., the 
average number of absorbed or emitted photons. For 
p« 1, only one-photon processes play a practical role. 
Within the single-sideband approximation AT is propor- 
tional to p, so that for small p the net transmission AT/p 
is independent of driving strength \U n \. 

The second relevant parameter is the ratio Huj/El be- 
tween the photon energy and the ballistic Thouless en- 
ergy El = hvp/L, i.e., the energy scale associated with 
the length L of the driven region, with vp the Fermi en- 
ergy. This scale governs the formation of Fabry-Perot 
scattering resonances (different from the polarization- 
dependent resonance to be identified below), so that 
Hlj/El provides a measure for the number of resonances 
within a window Hw. For Huj/El ^> 1, evanescent 
modes are densely spaced, which means that the quan- 
tum pumping mechanism by these modes is maximally 
efficient. 17 

We concentrate on the regime of universal and efficient 
quantum pumping of W ^> L, p <C 1 and El ^> Huj. The 
pumped current (3) then becomes independent of L, and 
the response to the driving is proportional to the effective 
driving strength p, Eq. (4). Extracting this factor and 
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FIG. 2. Total current driven across a wide and short graphene 
pump, as a function of Fermi energy (measured from the Dirac 
point in the static sheet) for various driving mechanisms as 
indicated in the legend. 



the scales of the energy and momentum integrals in the 
current formula, we obtain a characteristic current 
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A laser in the visible red, with wavelength of ~ 850 nm, 
1 mW output power and 1 fim spot size gives a typical 
coupling energy \U n \ = 0.1 me V, and a photon energy 
Huj = 0.7 eV. Hence p « 10~ 9 , and Huj/E l « 2200, 
such that the conditions for efficient pumping are ful- 
filled. The characteristic current is then Iq ~ 0.7 nA. For 
a far infrared laser with a wavelength of 50 /im, a power 
of 1 mW and a spot radius of 50 /im, we have p « 7 • 10 -6 
and Huj/El ~ 1900, and one may achieve much larger 
pumped currents of the order io ~ 50 nA. 

The precise value of the pumped current can be ob- 
tained by wave matching of modes with fixed transverse 
momentum k y and subsequent integration of the contri- 
butions of these modes. The resulting dependence on 
the Fermi energy (measured from the Dirac point in the 
static sheet) is shown in Fig. 2. As a consequence of 
particle-hole symmetry in the Dirac Hamiltonian, the 
pumped current is antisymmetry and vanishes at the 
Dirac point. Far from the Dirac point, i.e., for strongly 
doped graphene with \E F \ > Huj, the current saturates 
at a value of the order Io. Thus, Eq. (5) indeed gen- 
erally reflects the magnitude of the achievable pumped 
current. For a laser polarized in transport direction x, 
however, the pumped current displays a resonant fea- 
ture at \Ep\ = hu/2, where I ^> Iq. The singularity is 
logarithmic and grows as the ratio Huj/El is increased. 
Circular polarization (not shown) produces a very similar 
resonant response. In contrast, for transverse polariza- 
tion (along y), we witness only a small hump, while driv- 
ing by an ac gate voltage does not display any resonant 
features. 17 

This enhancement of the pumped current for x po- 
larization arises from a resonant coupling between the 
pseudo-spin states at Ep = ±fkj/2. This is revealed 
by inspecting the transmission imbalance AT(e, k y ), pre- 
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FIG. 3. Normalized net transmission AT/p, with p < 1 and 
Jtw/El = 50, for laser-driven graphene with polarization in x 
(longitudinal) and y (transverse) directions. The transverse 
momentum is measured in units of kuj = Huo/vf- The color 
scale (top bar) corresponds to positive current (left to right) 
in red, and negative in blue. 



sented in Fig. 3. A substantial pumping response is ob- 
tained in a diamond-shaped area, whose borders are de- 
fined by the transverse momenta at which modes in the 
sidebands n = 0, ±1 turn from propagating to evanes- 
cent. This delineates the evanescent-mode pumping 
mechanism which also dominates the response to an 
ac gate voltage. 12 ' 17 Notably, however, for laser driv- 
ing a strong pumping response also arises from the sub- 
harmonic resonances at energies e ~ These res- 
onances extend into the range of propagating incoming 
electrons with hvpk y < e. In the case of x-polarization 
the resonance evolves into a singular peak around the 
tips of the evanescent diamond, at vpk y = uo/2. 

This additional feature in the pumping response, which 
is absent for driving by ac gate voltages, reveals a fun- 
damental feature of the coupling of laser fields to the 
chiral charge carriers in graphene. For fixed wavevector 
fc, the static Hamiltonian in the pump, Hq = hvpk • tr, 
possesses propagating eigenstates |fc, s) with band index 
s = ±1, energies = shvp\k\ and pseudospin par- 
allel to sk. A laser with frequency uj induces resonant 
transitions between any two such states with the same 
k y , and energy difference \e kxkyjS - ek' x k y ,s>\ = ^j- The 
transition probability P = | (k' x , k y , s' | a ■ n | k x , k y , s) | 2 de- 
pends on the matrix element of the normalized driving 
operator U n /\U n \ = cr n. This probability becomes 
maximal, P = 1, if fc, k' JL n, and sk — —s'k. This res- 
onant condition can be satisfied for n in the x-direction, 
when vpk = (0,u;/2) and s = — s', which is manifested 
as a large pumping response at |e| ~ hvpk y « hu/2 (see 
Fig. 3, left panel). The response around this point has 
a constant positive (negative) sign for |e| below (above) 
Hlj/2 and thus builds up to a large peak in the total cur- 



rent. In contrast, for y polarization the response around 
e = ±hu/2 oscillates as a function of k y: which leads 
to a cancellation of different modes and a small overall 
contribution to the total pumped current. 

In summary, driving graphene by laser irradiation can 
lead to a sizable pumped current, which furthermore is 
resonantly enhanced when the polarization is parallel to 
the transport direction. For lasers in the visible range, 
typical currents are Iq = 0.7 nA, while the resonant con- 
dition \Ep\ = hu/2 corresponds to a relatively high dop- 
ing, which may require chemical functionalization or liq- 
uid electrodes. Moreover, infrared lasers may be em- 
ployed to lower the resonant condition to a much more 
accessible range of dopings (\Ef\ ~ 12meV), while at 
the same time increasing Iq and suppressing possible in- 
elastic effects due to electron-phonon interactions. The 
constraint in this case lies in keeping the laser spot size 
L below the mean free path in the flake. 
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